A two-layer model for studying 2D dissolved pollutant runoff over impermeable surfaces

Dissolved pollutants in stormwater are a main contributor to water pollution in urban environments. However, many existing transport models are semi-empirical and only consider one-dimensional flows, which limit their predictive capacity. Combining the shallow water and the advection–diffusion equations, a two-dimensional physically based model is developed for dissolved pollutant transport by adopting the concept of a ‘control layer’. A series of laboratory experiments has been conducted to validate the proposed model, taking into account the effects of buildings and intermittent rainfalls. The predictions are found to be in good agreement with experimental observations, which supports the assumption that the depth of the control layer is constant. Based on the validated model, a parametric study is conducted, focusing on the characteristics of the pollutant distribution and transport rate over the depth. The hyetograph, including the intensity, duration and intermittency, of rainfall event has a significant influence on the pollutant transport rates. The depth of the control layer, rainfall intensity, surface roughness and area length are dominant factors that affect the dissolved pollutant transport. Finally, several perspectives of the new pollutant transport model are discussed. This study contributes to an in-depth understanding of the dissolved pollutant transport processes on impermeable surfaces and urban stormwater management.     

Early solar system aqueous activity: K isotope evidence from Allende

The alkali element K is moderately volatile and fluid mobile; thus, it can be influenced by both primary processes (evaporation and recondensation) in the solar nebula and secondary processes (thermal and aqueous alteration) in the parent body. Since these primary and secondary processes would induce different isotopic fractionations, K isotopes could become a potential tracer to distinguish them. Using recently developed methods with improved precision (0.05‰, 95% confidence interval), we systematically measured the K isotopic compositions and major/trace elemental compositions of chondritic components (18 chondrules, 3 CAIs, 2 matrices, and 5 bulks) in the carbonaceous chondrite fall Allende. Among all the components analyzed in this study, CAIs, which formed initially under high‐temperature conditions in the solar nebula and were dominated by nominally K‐free refractory minerals, have the highest K₂O content (average 0.53 wt%) and have K isotope compositions most enriched in heavy isotopes (δ⁴¹K: ?0.30 to ?0.25‰). Such an observation is consistent with previous petrologic studies that show CAIs in Allende have undergone alkali enrichment during metasomatism. In contrast, chondrules contain lower K₂O content (0.003–0.17 wt%) and generally lighter K isotope compositions (δ⁴¹K: ?0.87‰ to ?0.24‰). The matrix and bulks are nearly identical in K₂O content and K isotope compositions (0.02–0.05 wt%; δ⁴¹K: ?0.62 to ? 0.46‰), which are, as expected, right in the middle of CAIs and chondrules. This strongly indicates that most of the chondritic components of Allende suffered aqueous alteration and their K isotopic compositions are the ramification of Allende parent‐body processing instead of primary nebular signatures. Nevertheless, we propose the small K isotope fractionations observed (< 1‰) among Allende components are likely similar to the overall range of K isotopic fractionation that occurred in nebular environment. Furthermore, the K isotope compositions seen in the components of Allende in this study are consistent with MC‐ICP‐MS analyses of the components in ordinary chondrites, which also show an absence of large (10‰) isotope fractionations. This is not expected as evaporation experiments in nebular conditions suggest there should be large K isotopic fractionations. Nevertheless, possible nebular processes such as chondrules back exchanging with ambient gas when they formed could explain this lack of large K isotopic variation.     

A nearly analytic symplectic partitioned Runge–Kutta method based on a locally one-dimensional technique for solving two-dimensional acoustic wave equations

In this paper, we develop a new nearly analytic symplectic partitioned Runge–Kutta method based on locally one-dimensional technique for numerically solving two-dimensional acoustic wave equations. We first split two-dimensional acoustic wave equation into the local one-dimensional equations and transform each of the split equations into a Hamiltonian system. Then, we use both a nearly analytic discrete operator and a central difference operator to approximate the high-order spatial differential operators, which implies the symmetry of the discretized spatial differential operators, and we employ the partitioned second-order symplectic Runge–Kutta method to numerically solve the resulted semi-discrete Hamiltonian ordinary differential equations, which results in fully discretized scheme is symplectic unlike conventional nearly analytic symplectic partitioned Runge–Kutta methods. Theoretical analyses show that the nearly analytic symplectic partitioned Runge–Kutta method based on locally one-dimensional technique exhibits great higher stability limits and less numerical dispersion than the nearly analytic symplectic partitioned Runge–Kutta method. Numerical experiments are conducted to verify advantages of the nearly analytic symplectic partitioned Runge–Kutta method based on locally one-dimensional technique, such as their computational efficiency, stability, numerical dispersion and long-term calculation capability.     

结合变差纹理特征的极化SAR建筑物震害信息提取

快速评估建筑物地震灾害信息对地震应急救援工作有着指导意义,而极化SAR具有全天候、全天时的特点,因此利用极化SAR图像提取震害信息已逐渐成为研究热点。虽然极化SAR具有丰富的极化信息,但其纹理信息不可忽略,尤其是完好的人工建筑物在图像上呈现规则的纹理特征,而倒塌建筑区域纹理分布杂乱,因此结合纹理信息也可以很好地提取建筑物信息。以2010年玉树地区的全极化SAR数据为研究对象,首先,利用Yamaguchi分解的体散射分量PV提取了SAR图像中的建筑物区域以及道路、水系等非建筑物信息,在此基础上,对相干散射矩阵T11分量中倒塌建筑物、完好建筑区域进行变差计算,根据变差曲线确定变程a后,再对建筑物区域采取窗口m*m(m=3*a)进行变差计算得到变差纹理信息,最后利用FCM算法对变差纹理信息分别提取完好建筑物和倒塌建筑物区域,为了对比分析,文章利用Yamaguchi分解的二次散射分量PD提取完好建筑物区域,与震后光学遥感图像对应样本点进行人工验证,得到完好建筑物的提取精度为80.18%,倒塌建筑物的提取精度为84.54%,道路水系的提取精度为77.58%。     

高分六号宽幅相机在轨几何定标及精度验证

高分六号宽幅相机能够实现单相机成像幅宽优于800 km,对大尺度地表观测和环境监测具有独特优势。在轨几何定标是光学遥感卫星几何处理的关键环节,直接影响影像的几何质量。本文充分考虑高分六号宽幅相机超大视场的畸变特性以及多谱段的成像特点,提出宽幅相机在轨几何定标方法,采用基于探元指向角的几何定标模型补偿宽幅相机系统误差,通过绝对定标和相对定标方法联合估计各波段的内外定标参数。利用Landsat 8影像、资源三号DSM为参考数据,对宽幅相机进行绝对定标处理,再利用ASTER GDEM为参考数据进行相对定标处理,其几何定标结果表明,高分六号宽幅影像绝对定位精度在3像素左右,内部几何精度能稳定在1像素,且波段间配准精度在0.3像素以内,表明在轨几何定标后高分六号宽幅影像几何质量得到了明显的提升。     

滨海城市竖向避难场所选址可靠性评价与若干问题探讨

关注滨海城市竖向避难场所的选址可靠性问题，首先深化确立“竖向避难场所”概念；然后构建竖向避难场所选址可靠性评价模型，综合自然地理要素和建成环境要素，从与灾害直接作用相关的选址暴露性、与应急交通相关的选址敏感性、与应急服务相关的选址适应性3个维度展开，搭建3级评价指标体系；再以天津滨海新区为例、聚焦滨海城市潮洪灾害和人口安全矛盾的集核——临港城区展开实证研究，借助ArcGIS分类与可视化发现：高低可靠性选址之间具有明显的空间分异，可以分解为灾害暴露性的“近岸高?远岸低”分异、交通敏感性的“外围高?中心低”分异、服务适应性的“中心与沿河高?外围低”兼“近港高?远港低”分异。同时发现：商业设施用地选址价值较高，在高可靠性选址中占比第一；中小学和社会福利设施用地选址价值最高，在高可靠性选址中占比第二；文化科研和娱乐康体设施用地在高可靠性选址中占比最小。针对临港城区实证研究结果，提出滨海城市竖向避难场所选址及可靠性提升对策。     

生物土壤结皮对风沙土和黄绵土膨胀特性的影响

作为重要的土壤物理性质，膨胀性在影响土壤导水性、持水性、抗蚀性以及土壤结构的形成和发育等方面发挥着重要作用。为了探讨生物土壤结皮（BSCs）土壤的膨胀特性及其主要影响因素，针对黄土高原风沙土和黄绵土两种典型土壤，利用膨胀仪测定并比较了有、无藓结皮及其在不同因素（初始含水量、干湿循环、冻融循环、温度）下膨胀率的差异，分析了BSCs对土壤膨胀性的影响及其与环境因素和BSCs性质的关系。结果显示：风沙土上藓结皮的膨胀率为1.93%，较无结皮增加了8.65倍；而黄绵土上藓结皮的膨胀率为2.05%，与无结皮相比降低了76.68%。藓结皮的生物量和厚度与其膨胀率在风沙土上均呈线性正相关关系（P < 0.05），在黄绵土上分别呈二次函数（P=0.02）和线性正相关关系（P=0.02）。初始含水量同时影响了土壤最大膨胀率和稳定膨胀时间，影响程度风沙土远大于黄绵土（包括藓结皮和无结皮）；干湿循环次数对无结皮土壤膨胀率的影响程度大于藓结皮土壤，其中风沙土和黄绵土上无结皮的膨胀率分别是50.00%~620.00%和-2.28%~10.81%，而两种土壤上藓结皮的膨胀率分别是-5.70%~10.88%和-10.24%~-21.46%；冻融循环下4种土壤的膨胀率均有不同程度的降低，降幅为0~18.54%。黄绵土无结皮的膨胀率受温度影响程度较大，50℃下黄绵土无结皮的膨胀率分别是25℃和35℃下的1.17倍和1.21倍。BSCs显著地改变了风沙土和黄绵土表层的膨胀性，其影响的程度和方向取决于土壤类型。同时，BSCs的膨胀性受含水量、温度、干湿以及冻融循环等关键因素影响。     

游客-环境互动视角下旅游者目的地环境恢复性效应

日常生活环境和旅游目的地环境都可以为个体提供恢复的机会，不同环境的属性特征导致不同程度的恢复。为了探究环境类型与恢复性效应之间的关系及停留时间的调节作用，文章采用问卷调查方式，利用配对样本t检验和混合设计方差分析相结合的方法，探讨游客在不同环境中的恢复性效应，并对这一现象提供心理上的解释。结果表明：1）恢复体验的4个维度（放松体验、心理脱离、掌握体验、控制体验）在旅游目的地环境中的得分均高于日常生活环境。2）掌握和控制体验维度在城市旅游目的地中的变化程度高于乡村旅游目的地，两种类型目的地的放松和心理脱离维度差异较小。3）随着在目的地停留时间的延长，乡村旅游目的地放松体验和心理脱离维度变化程度高于城市旅游目的地，两种类型目的地环境恢复性效应会相互靠近。因此，游客的环境恢复性效应不仅受到环境和目的地类型的影响，同时还与游客的适应程度有关，且受到其停留时间的调节作用。     

Salinivibrio sp.YH4胞外丝氨酸蛋白酶EYHS耐盐性及生物信息学分析

摘要：目的 探讨菌株Salinivibrio sp.YH4分泌的丝氨酸蛋白酶EYHS的耐盐性及结构特征。方法 明胶底物酶谱法分析EYHS的耐盐性。应用生物信息学手段对EYHS及6种耐盐的S8家族丝氨酸蛋白酶结构特征进行分析。结果 EYHS在4 mol/L的NaCl溶液中仍具有活性，属于耐盐蛋白酶。EYHS及6种S8家族丝氨酸蛋白酶分子表面的loop区等无规则卷曲所占比例较高，α-螺旋与β-片层则主要位于酶分子内部。EYHS分子表面酸性氨基酸含量较高，且具有弱疏水内核。多序列比对发现蛋白酶的催化三联体两侧存在高度保守的基序和保守的极性氨基酸及芳香族氨基酸，并存在多个保守的Gly与Ala。同源模建和表面电荷分布显示，α螺旋和β片层围成了蛋白酶的催化腔，EYHS活性中心包含由Asp32、His65与Ser215组成的催化三联体，且催化位点区域表面静电势为负。结论 上述结构特征可能有助于耐盐丝氨酸蛋白酶EYHS在高盐环境下维持其稳定性和适度柔性，并有助于其催化功能的发挥，为深入研究耐盐丝氨酸蛋白酶的高盐环境适应性提供了一定的理论依据。     

3D Mineral Prospectivity Mapping with Random Forests: A Case Study of Tongling,Anhui, China

Natural Resources Research - In the past few decades, a variety of data-driven predictive modeling techniques has led to a dramatic advancement in mineral prospectivity mapping (MPM). The random...